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Abstract. The use of radioactive ion beams is shown to offer the possibility to study collective pairing 
states at high excitation energy, which are not usually accessible with stable projectiles because of large 
energy mismatch. In the case of two-neutron stripping reactions induced by ^He, we predict a population 
of the Giant Pairing Vibration in ^°^Pb or '^^^Sn with cross sections of the order of a millibarn, dominating 
over the mismatched transition to the ground state. 

PACS. 21.60.Ev Collective Models - 25.60. Je Transfer Reactions 



1 Introduction. 

Large efforts have been recently dedicated to the study 
of different aspects of reaction mechanism in collisions 
induced by weakly-bound radioactive beams. The long 
tails of the one-particle transfer form factors due to the 
weak binding, associated with the possibility of unusual 
' behaviour of pairing interaction in diluted systems, has 
raised novel interest in the possibility of studying the pair 
field via two-particle transfer processes with unstable beams 
]l|. On the other hand, in transfer reactions induced by 
■ weakly bound projectiles on stable targets, the Q-values 
' for the low-lying states will be very large (typically of 
the order of 10-15 MeV for the (^He,'*He) stripping re- 
action). This will strongly hinder these processes for re- 
actions where the semi-classical optimum matching condi- 
tions apply, as it is the case of bombarding energies around 
the Coulomb barrier on heavy target nuclei. Higher bom- 
barding energies, where the matching conditions are less 
stringent, may on the other hand not be suitable because 
of large break-up cross sections. The same matching con- 
ditions will favour instead the population of highly ex- 
cited states, as the Giant Pairing Vibrations (GPV), and 
the use of Radioactive Ion Beams (RIB) may therefore 
become instrumental in offering the opportunity of study- 
ing nuclear structure aspects that are not usually acces- 
sible with stable projectiles. These Giant Pairing Vibra- 
tions are in fact predicted to have strong collective 
features, but their observation may have so far failed 
because of large mismatch in reactions induced by pro- 
tons or tritons, at variance to the case of the low-lying 
pairing vibrations, which have been intensively and suc- 
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cessfuUy studied around closed shell nuclei in two-particle 
transfer reactions [Q. All these 0+ states are associated 
with vibrations of the Fermi surface and are described in 
a microscopic basis of the shell model as correlated two 
particle- two hole states. In the case of the Giant Pair- 
ing Vibrations the excitation involves the promotion of a 
pair of particles (or holes) in the next major shell (hence 
an excitation energy around 2fiuj) and is expected to dis- 
play a collective pairing strength comparable with the low- 
lying vibrations. Also in the case of superfluid systems in 
an open shell the system is expected to display a collec- 
tive high-lying state, that in this case collects its strength 
from the unperturbed two-quasiparticle 0"*' states with en- 
ergy 2hLu. To investigate this possibility we made estimates 
of cross sections to the Giant Pairing Vibrations in two- 
particle transfer reactions, comparing the cases of bound 
or weakly-bound projectiles. As examples we have consid- 
ered the case of (^''C,^^C), from one side, and the case 
of (^IIe,''He) as representative of a reaction induced by a 
weakly-bound ion. As targets, we have chosen the popu- 
lar cases of the lead and tin regions (so considering both 
"normal" and "superfluid" nuclei). To perform the calcu- 
lation, we will first evaluate the response to the pairing 
operator in the RPA, including both the low-lying and 
high-lying pairing vibrations. As a following step we will 
then construct two-neutron transfer form factors, using 
the "macroscopic" model for pair-transfer processes. Fi- 
nally, estimates of cross sections will be given using stan- 
dard DWBA techniques. As we will see, in the case of the 
stripping reaction induced by ^He, the population of the 
GPV is expected to display cross sections of the order of a 
millibarn, dominating over the mismatched transition to 
the ground state. 
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The paper is organized as follows. In the next sec- 
tion we discuss the theoretical formalism used for normal 
and for superfluid nuclei. In section 3 we recall the basics 
aspects of the macroscopic form factors for two-particle 
transfer reactions and in section 4 we display the results 
of calculations for the paradigmatic examples of ^osp^ g^j^fj 
ii^Sn. 



2 The pairing response and the Giant Pairing 
Vibration. 

A simple way of displaying the amount of pairing corre- 
lations is in terms of the pair transfer transition densities 
1^. These are defined as the matrix element of the pair 
density operator connecting the ground state in nucleus 
A with the generic 0+ state |n) in nucleus A ±2, namely 

6pp{r) = {n\pp\0), (1) 

where the generalized density operator is given by 



VWTT 
An 



Here Ra{r) are the radial wave functions of the a — {nlj} 
level and the sum runs over both particle and hole lev- 
els. The pair transfer strength to each final state can be 
obtained from the corresponding pair transfer transition 
density by simple quadrature, namely 



f3p — / Anr^Sppdr. 



(3) 



For normal systems around closed shell the strong L=0 
transition follows a vibrational scheme, where the corre- 
lated pair of fermions (pairing phonon) change by one Q . 
In this case, there are two types of phonons associated 
with the stripping and pick-up reactions. The two-particle 
collective state is called "addition" pairing phonon while 
the two-holes correlated state is known as " removal" pair- 
ing phonon. From a microscopic point of view the two 
kind of phonons, corresponding to the {A ± 2) nuclei can 
be described in terms of the two-particle (two-hole) states 
of the Tamm DancofF Approximation(TDA) or in a better 
way by a Random Phase Approximation (RPA) . We start 
from a hamiltonian with a Monopole Pairing interaction: 



where 
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Here the Oj creates a particle in an orbital j, where j 
stands for all the needed quantum numbers of the level. 
The constant G is the strength of the pairing interaction 
and the coefficients M{ji,j2) are: 



M{ji,.n) 
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(6) 



where the detailed radial dependence of f{r) is taken to 
be of the form and in our case is a constant since we 
are dealing only with L = states. The Pairing phonons 
are defined for closed shell nuclei as: 
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where k{i) stands for levels above (below) the Fermi level. 
We will use j for any of both cases. X„ and y„ are the 
forward and backward amplitudes. Within this model the 
pair transfer strength associated with each RPA state is 
microscopically given by 



(8) 



In Fig. la we display the predicted pairing response in 
the case of ^°^Pb, namely two-neutron holes with respect 
to the double magic ^°*Pb. The set of single-particle levels 
that has been used in the RPA calculation, was obtained 
using the spherical harmonic oscillator levels with correc- 
tions due to the centrifugal and spin-orbit interactions 
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where htu = 41A~3, A is the mass number of the nu- 
cleus, N is the principal quantum number and j, I are the 
total and orbital angular momentum quantum numbers, 
respectively. The quantities k and p are parameters cho- 
sen to obtain the best fit for each nucleus ||]. We have 
included in the calculation all the single-particle levels 
starting from = up to 10. This set is expected to 
be good enough for our calculation of the Giant Pairing 
Resonance, except for the levels around the Fermi surface. 
In the lead region we prefer to use experimental values for 
the shells just above and below the Fermi surface |jl^ . The 
Figure shows, in addition to the strong collectivity associ- 
ated with the ground state transition, a strong collective 
state with about half of the g.s. strength at high excita- 
tion energy, around 16 MeV, which can be interpreted as 
the Giant Pairing Vibration. Similar situation is shown in 
Fig. lb for the corresponding two-neutron addition states 
in the ^^"^Pb. Again one may interpret the strength at 
about 12 MeV as associated with the giant mode. 

We consider now the case of superfluid spherical- nuclei. 
In this case we make a BCS transformation of the hamilto- 
nian defined in Eq. Q changing from particle to quasipar- 
ticle operators, introducing the usual occupation param- 
eters. We start from a single-quasiparticle Hamiltonian 
plus a two-quasiparticle interaction corresponding to the 
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Fig. 1. Pairing response for the removal (upper panel) and 
addition (lower panel) mode in •^"**Pb. The ground state tran- 
sition and the candidate for GPV are evidenced. 



residual H22 + H^q of the pairing force 

3 

+ 27rG^M(ii,ii)M(j2,j2)- 



where 



Uf = 



(10) 

(11) 
(12) 



v. 



(13) 



The energies Ej = ^e'^ + are the quasi-particle ener- 
gies, and e j = e — A arc the singlc-particlc energies with 
respect to the chemical potential A and A is the BCS gap. 



3-m, 



As usual we have defined aj = aj:^ = (—1) 

For superfluid systems the addition and removal RPA 
phonons cannot be treated separately. The dispersion re- 
lation, that relates the strength of the interaction with the 
energy-roots of the RPA, becomes a two by two determi- 
nant. Prom the RPA equations: 

rl = Y,{^n{j)[a}a%+Ynij)[ajaj]oo) (14) 



we can obtain the following factors 
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and the dispersion relation is in this case: 
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It can be shown that a; = is solution of that equation 

and correspond to the Goldstone boson corresponding to 
the breaking of the number of particle symmetry. Once we 
have obtained the energies Un of the different RPA roots, 
we can write the components of the RPA phonon in the 
form: 



Xn{j,j) = 



47rGM(j,j) 



Ej + Ej 

4^GM(j,j 



uf 



A-kGz 



Ej +Ej+uJr,\ ^ (1 - 47rGy) 



+ ^/ ) ^n, 



(20) 



where yl„ is determined by normalizing the phonon corre- 
sponding to the n— th root of the RPA. The normalization 
condition reads 



E[^^(j)-^„'(i)] = i- 



(21) 
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Finally, we can obtain for each state n the pairing strength 
parameter (3p with the following formulae: 

/3p(2p) = V2jTT(n|[a]a]]oo|0) = 

3 

3 

(3p{2h) ^ J2 V27TT(n|[a,a,]oo|0) = 

3 

= ^2jTl[VfX^{j) - U^Y^m (22) 



The predictions of the pairing strength distribution for the 
superfluid system ^^^Sn are shown in the two panels of Fig. 
2. For the calculation we have used the single-particle lev- 
els from Ref. |l§ . These last ones have been proved to give 
good results in BCS calculations using a pairing strength 
G = g/A, where g ~ 20MeV. We assume that the rest of 
the levels have occupation probability 1(0) if they are far 
below(above) the Fermi surface. The change of the single 
particle energies around the Fermi surface has been done, 
in both cases, taking care of keeping the energy-centroids 
of the exchanged levels in the same position. The figure 
clearly shows the occurence of high-lying strength which 
can be associated with the Giant Pairing Vibration. Note 
that, with respect to the case of ^°*Pb, there is a minor 
fragmentation of the strength both in the low-lying and 
in the high-lying energy region. 



3 Macroscopic form factors for two-particle 
transfer reactions. 

The description of the reaction mechanism associated with 
the transfer of a pair of particles in heavy ion reactions 
has always been a rather complex issue. In the limit in 
which the field responsible for the transfer process is the 
one-body field generated by one of the partners of the re- 
actions, at least for simple configurations the leading order 
process is the successive transfer of single particles. In this 
framework the collective features induced by the pairing 
interaction arise from the coherence of different paths in 
the intermediate (A+1 , A-1) channel due to the corre- 
lation present in the final (A+2) and (A-2) states. The 
actual implementation of such a scheme may turn out not 
to be a simple task, due to the large number of active 
intermediate states, and the use of a simpler approach 
may be desirable. This is offered, for example, by the 
"macroscopic model" for two-particle transfer reactions, 
that parallels the formalism used to describe the inelastic 
excitation of collective surface modes. In that case, as an 
alternative to the (more correct) microscopic description 
based on a superposition of particle-hole excitations, one 
has traditionally resorted to collective form factors of the 
form g 

Fx{r) = PxR^, (23) 
dr 
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Fig. 2. Pairing response for the removal (upper panel) and ad- 
dition (lower panel) mode in ^^^Sn. The ground state transition 
and the candidate for GPV are evidenced. 



in terms of the radial variation of the ion-ion optical poten- 
tial U induced by the surface vibrations, with the strength 
parameter f}\ obtained from the strength of the B(EA) 
transition. In the case of the pair transfer, the correspond- 
ing vibration is the fluctuation of the Fermi surface with 
respect to the change in the number A of particles, and 
the corresponding form factor Fp is assumed to have the 
parallel form [|j 



Fp{r) 



dU 



(24) 



in terms of the "pairing deformation" parameter j3p asso- 
ciated with that particular transition, defined in the previ- 
ous section. The assumption of simple scaling law between 
nuclear radius R and mass number A allows to rewrite the 
two-particle transfer form factor into an expression which 
is formally equivalent to the one for inelastic excitation, 
namely 
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This formalism has been successfuUy applied to quite a 
number of two-particle transfer reactions . As in the 

case of inelastic excitations, macroscopic collective form- 
factors may in some cases only give a rough estimate to the 
data, requiring more elaborate microscopic descriptions. 
Nonetheless, the use of simple macroscopic form factors 
is of unquestionable usefulness in making predictions, in 
particular in cases, as the one we are discussing, where 
experimental data are not yet available and estimates are 
needed in order to plan future experiments. 



4 Applications: estimates of two-neutron 
transfer cross sections. 

In order to evidence the possible role of unstable beams in 
the study of high-lying pairing states, we compare in this 
section two-particle transfer reactions induced either by a 
traditionally available beam (e.g. the (^''C, ^^C)) or by a 
more exotic beam (e.g. the reaction (^He, '^He)). As a tar- 
get, we have considered the two cases of ^°^Pb and ^^^Sn, 
as representative cases of normal and superfluid systems 
in the pairing channels. In both cases, we have consid- 
ered the full pairing L=0 response, e.g. all transitions to 
0+ states in ^lOpb and "^Sn, as described in Sect 2. The 
Q-values corresponding to the transitions to the ground- 
states and to the GPV states are displayed in Table 1. 
For each considered state the two-particle transfer cross 



a) 







''He He 


^^'^Sn- 


^ ^^«Sn,, 


3.15 MeV 


15.298 MeV 


208pb _ 




-4 MeV 


8.148 MeV 


"«Sn^ 


^""^^^SncFv 


-6.746 MeV 


5.402 MeV 


208 ^ 


fDGPV 


-15.81 MeV 


-3.662 MeV 



Table 1. Q- values for ground-state and GPV transitions. The 
target (column) and projectile (row) are specified. 



section has been calculated on the basis of the DWBA 
(using the code Ptolemy Q) employing the macroscopic 
form factor described above, with a strength parameter 
as resulting from the RPA calculation. For the ion-ion op- 
tical potential, the standard parameterization of Akyuz- 
Winther |l^ has been used for the real part, with an imag- 
inary part with the same geometry and half its strength. 
In all cases, the bombarding energy has been chosen in 
order to correspond, in the center of mass frame, to about 
50% over the Coulomb barrier. The angle-integrated L=0 
excitation function is shown in Fig. 3b as a function of 
the excitation energy for the 208pb(i4c,i2c)208pb re- 
action at Ecm=95 MeV. For a more realistic display of 
the results, the contribution of each discrete RPA state is 
distributed over a lorentzian with r= k E^, with k ad- 
justed to yield a width of AMeV for the giant pairing vi- 
bration. As the figure shows, the large (negative) Q-value 
associated with the region of the GPV (see Table 1) com- 
pletely damps its contribution, and the excitation function 
is completely dominated by the transition to the ground 
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Fig. 3. Differential cross-sections as a function of the exci- 
tation energy for the two reactions : a) ^°**Pb(''He,''He)^°*Pb, 
and b) 2'^**Pb("C,'^C)^'^**Pb. See text for details. 



state and the other low-lying states. The situation is very 
different for the 208pb(6He,''He)208pb reaction at Ec„i=41 
MeV, whose excitation function is shown in Fig. 3a. In this 
case the weak binding nature of ^He projectile leads to a 
mismatched (positive) Q-value for the ground-state tran- 
sition {Qgs= 8.148 MeV), favouring the transfer process 
to the high-lying part of the pairing response. In this case 
the figure shows that, in spite of a smaller pairing matrix 
element, the transition to the GPV is of the same order 
of magnitude of the ground-state transfer (1.8 mb for g.s. 
and 3.1 mb for the GPV). Note that a total cross section 
to the GPV region of the order of some millibarn should 
be accessible with the new large-scale particle-gamma de- 
tection systems. 

A similar behaviour is obtained in the case of a tin tar- 
get. In Figs. 4a and 4b the corresponding excitation func- 
tions for the "6sn(i4c,i2c)ii8sn reaction (at Ec,„=69 
MeV) and the "6Sn(6He,''He)"SSn reaction (at Ec™=40 
MeV) are compared. Now the transition to the GPV dom- 
inates over the ground-state transition when using an He 
beam ( 0.4 mb for g.s. and 2.4 mb for the GPV). From 
a comparison with the RPA strength distributions of Fig. 
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Fig. 4. Differential cross-sections as a function of the excitar- 
tion energy for the two reactions : a) ^^^Sn(^He,*He)'^'^®Sn, and 
b) "«Sn("C,i2C)"*Sn. The comparison between the GPV 
and the ground-state clearly shows the different strength. No- 
tice the different vertical scale with respect to figure 3. 



1 and 2 one can see that the giant pairing vibrations is 

definitely favoured by the use of an ^He beam instead of 
the more conventional ^^C one, because the transition to 
the ground-state is hindered, while the GPV is enhanced 
(or not changed), because of the effect of the Q- value. 



5 Conclusions. 

The role of radioactive ion beams for studying different 
features of the pairing degree of freedom via two-particle 
transfer reactions is underlined. A ^He beam may allow an 
experimental study of high-lying collective pairing states, 
that have been theoretically predicted, but never seen 
in measured spectra, because of previously imfavourable 
matching conditions. The modification in the reaction Q- 
value, when passing from ^^C to ^He, that is a direct con- 
sequence of the weak-binding nature of the latter neutron- 
rich nucleus, is the reason of the enhancement of the tran- 
sition to the giant pairing vibration with respect to the 
ground-state. 



